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Abstract. We present in this article two different ways to make CMB maps in 
practice, from large timelines. One is to make a simple destriping, fitting the data 
and using the scan intercepts to remove the low frequency noise (stripes). The sec- 
ond, optimal, is to resolve linearly the map-making problem, which in case of big 
timelines must be simplified and changed from matrices to vectors for the calcula- 
tions. Assuming few conditions on the noise, it is possible to make fast map-making 
tools. 



1 Introduction 

The last ten and the future ten years provide a large number of CMB (Cos- 
mic Microwave Background) experiments. Their main goal is to estimate the 
cosmological parameters through sky temperature maps. The map-making is 
therefore a key question of these experiments, moreover with the advent of 
large time-ordered data (TOD) being currently analysed or simulated. The 
"brute force" direct inversion of the linear map-making problem is beyond 
the reach of computing facilities of today. Thus the data analysers need to 
find faster methods, which have to be as optimal in theory and as efficient in 
practice as possible. 

2 Destriping for Planck HFI data 

This is used to clean Planck High Frequency Instrument (Bersanelli et al. 
1996) simulated data, to be published in Giard et al. (2000). It uses a simple 
algorithm to destripe the data and make the map. The scanning strategy 
assumed to make the timeline is 1 round per minute of the beam on the 
sky, so that the timeline is a succession of circles on the sky. The cleaning 
algorithm follow the processing described below: First each circle is adjusted 
so that the measured signal fits a cosecant law (galactic dust), plus dipole 
emission: 

M,k = rfe {Dipole^k + — ^TTT + Cl) (1) 
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where k is the index of each circle of the data, i the data index along 
the circle, and the detector response. The free parameters are gk (galac- 
tic emission for the circle k) and C° (constant to be subtracted). Then we 
destripe the data using an algorithm derived from Delabrouille (1998) which 
uses the scan intercepts. We adjust the constants to subtract to each circle 
k by minimising the spread between the measurements contributing to the 
same sky pixel. 

3 The optimal method in practice 

3.1 The map-making problem 

If we assume the instrument to be linear, the time-ordered data, represented 
by a vector y, is linearly dependent on the real sky, represented by a vector 
X. We have thus: 

y ^ Ax + n (2) 

where A is the point-spread (convolution) matrix and n the random noise 
vector in the timeline. For CMB experiments, x would represent the pixelised 
sky map of the CMB temperature. The map-making problem is thus written 
as this: 

xo = Wy (3) 

where xq is the vector of the reconstructed sky map, and W the inversion 
matrix. The inversion matrix W depends on the method used. The simplest is 
the pixel averaging, but the optimal methods for estimating the map are the 
COBE method (with no prior for x) or the Wiener filter (with gaussian prior 
for x). See Tegmark (1997) for a set of linear and non-linear map-making 
methods. As an example, the COBE method is written as this: 

W ^[A'N-^A]-^A*N-^ (4) 
with N the noise covariance matrix N —< nn} >. 

3.2 The large timeline problem 

The CMB experiments now in analysis (Boomerang, MAXIMA...) or to come 
(Archeops, MAP, Planck...) provide several megabytes of time-ordered data. 
The noise covariance matrix, for instance, would be several terabytes of data, 
which is impossible to handle or even to write. How to solve optimally the 
map-making problem without writing any matrix anywhere ? 

Assuming no beam and the stationarity of the noise, it is possible to make 
the map iteratively handling only vectors. Using the COBE method (eq. 
the deconvolution equation (eq. 0) can be written as this: 
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[A'N~^A]xa = A'N-^y (5) 

If the noise is stationary, the matrix multiplication N^^ y is equivalent 
to a convolution by a kernel, i.e. a multiplication in the Fourier space. If the 
beam is not taken into account, A is making a timeline from a map, and A* is 
the pixel averaging of a timeline into a map. Eq. |^ leads to iterative methods 
such as: 

Xn+l ^ Xn +T1 A' N~^y~iq A' N^^AXn (6) 

where 77 is a free parameter of the iterative method. It is also possible to 
make converge the noise map instead of directly the sky map, by changing 
the variable xq into nx^ = xq — [A* A\~^ A*y. 




Fig. 1. Convergence for three pixels of Archeops technological flight data 

Fig. |l| shows the convergence of such methods, here applied to data from 
Archeops technological flight, realised at Trapani in sept. 99 (Benoit et al. 
2000). 
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